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ABSTRACT 

Context. Highly-collimated outflows are believed to be the earliest stage in outflow evolution, so their study is essential for understanding the 
processes driving outflows. The BHR71 Bok globule is known to harbour such a highly-collimated outflow, which is powered by a protostar 
belonging to a protobinary system. 

Aims. We aimed at investigating the interaction of collimated outflows with the ambient molecular cloud by using molecular tracers. 
Methods. We mapped the BHR71 highly-collimated outflow in CO(3-2) with the APEX telescope, and observed several bright points of the 
outflow in the molecular transitions CO(4-3), 13 CO(3-2), C I8 0(3-2), and CH30H(7-6). We use an LVG code to characterise the temperature 
enhancements in these regions. 

Results. In our CO(3-2) map, the second outflow driven by IRS2, which is the second source of the binary system, is completely revealed and 
shown to be bipolar. We also measure temperature enhancements in the lobes. The CO and methanol LVG modelling points to temperatures 
between 30 and 50 K in the IRS 1 outflow, while the IRS2 outflow seems to be warmer (up to 300 K). 
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1. Introduction 

Low-mass stars form in molecular clouds from the collapse of 
dense cores. In the earliest stage, known as Class 0, the newly- 
born protostar is deeply embedded in a thick envelope of gas 
and dust from which it will accrete most of its mass, and it 
drives powerful outflows. However, the mechanisms driving 
bipolar outflows, as well as the interaction of the outflow with 
the surrounding material, are still not understood well. It is 
commonly believed that highly-collimated outflows represent 
the earliest stage in outflow evolution, making their observa- 
tion central to getting new insights into the mechanisms driving 
outflows. 

The BHR71 isolated Bok globule, located at ~ 200 pc, har- 
bours one of the best examples of such a highly-collimated out- 
flow. This outflow, which is almost lying in the plane of the sky, 
wa s first discove r ed and mapped in I2 CO(1-0) and I3 CO(1-0) 
by iBourke et alJ Jl997l) . It appears to be driven by the very 
young stellar o bject IRAS 1 1590-6452, classified as a Class 
0/1 protostar by Froebrich (2005 ; j), with a total luminosity of 
~ 9 L . The CO emission was observed to peak towards the 
outflow lobes, implying that the ambient gas is heated by inter- 
action with the outflow. 

Subsequent observations by iGarav et"afl Jl998l) aimed at 
constraining the magnitude of the temperature enhancement in 



the lobes, by mapping several transitions of CS, SiO, CH3OH, 
and HCO + with the SEST telescope. Assuming n= 10 5 cm -3 , 
the CH3OH (2-1) and (3-2) observations were found to be con- 
sistent with kinetic temperatures in the 50- 1 00 K range. 

ISOCAM observations of BH R71 showed that the IRAS 
source is in fact a multiple system dMvers & Mardonesll 998). 
Two protostellar objects IRS 1 and IRS2 are present in BHR7 1 , 
and t he large collimated outflow is driven by IRS1 (Bourke 
2001). Combination of CO(2-l), NIR, and cm observations 
shows evidence of the presenc e of the blue l obe of a second, 
fainter outflow driven by IRS2 dBourkel200ll) . 

In this Letter, we present CO(3-2) observations that fully 
unveil the presence of this second outflow, and characterise the 
temperature enhancement towards bright knots of the larger 
outflow by observations of CO(4-3), 13 CO(3-2), 18 CO(3-2), 
and CH 3 OH(7-6). 



2. Observations 

Observations towards BHR71 were carried out in September 
and November 2005 with the Atacama Pathfinder Experiment 
(APEX 1 , Giisten et al. 2006, this volume), using the facility's 
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1 This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX). APEX is a collaboration between 
the Max-Planck-Institut fur Radioastronomie, the European Southern 
Observatory, and the Onsala Space Observatory. 
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Fig. 1. (a) CO(3-2) channel maps (centred on a(2000) = 12 h 01 m 37 s , (5(2000) = -65°08'53.5") ona/T; dv scale (the v lsl - of the 
cloud is -4.5 km s -1 ). (b) Map of the integrated CO(3-2) emission in the -14<v<-6 (blue) and -3<v<5 kms -1 (red) intervals, 
enlightening the two bipolar outflows. Contours are from 7.5 to 30 in steps of 4 Kkms -1 . Contours for higher flux (in the 
north-south outflow) have not been plotted in order to show the fainter outflow better. The two black stars show the position of 
the two protostellar sources. The black circle in the lower right corner shows the beam size, (c) CO(3-2), 13 CO(3-2), C 18 0(3-2), 
and CO(4-3) observations on the three positions marked as black squares in (b). For clarity, CO(3-2) and CO(4-3) lines have been 
shifted by 8 and 4, respectively, along the y-axis. 



own APEX2a receiver (Risacher et al. 2006, this volume) and 
the MPIfR First Light APEX Submillimeter Heterodyne re- 
ceiver (FLASH, Heyminck et al. 2006, this volume), along 
with the MPIfR Fast Fourier Transform Spectrometer backend 
(FFTS, Klein et al. 2006, this volume). The beam of the tele- 
scope is 17.3" at 345GHz and 13.3" at 460GHz (Giisten et 
al. 2006, this volume). Focus was checked at the beginning of 
each run and after sunrise on Saturn. Pointing was checked on 
Saturn, 07454-7112 or IRAS 15 194-51 15. The pointing accu- 
racy was found to be of the order of 4" rms. 



We mapped the large-scale outflow in the CO(3-2) transi- 
tion using the on-the-fly observing mode. Several other transi- 
tions were observed towards several bright points of the out- 
flow (CO(4-3), 13 CO(3-2), C ls O(3-2), and the CH 3 OH 7* ^6 k 
band). The double-side-band system temperatures ranged from 
100 K to 200 K for all observations, except for CO(4-3) where 
they were in the range 1 000-3000 K. We used 8192 channels of 
the FFTS over the 1 GHz bandwidth, yielding an effective res- 
olution of about 240 kHz. The observations were performed in 
position-switching mode using the APECS software (Muders 
et al. 2006, this volume). The data were reduced with the 
CLASS software (see http://www.iram.fr/IRAMFR/GILDASl. 
Antenna temperatures were converted to main-beam tempera- 
tures using forward and beam efficiencies, respectively, of 0.97 
and 0.73 at 345 GHz and of 0.95 and 0.61 at 464 GHz (Glisten 
et al. 2006, this volume). 



3. Analysis 

3.1. Outflow morphology 

The channel map of the CO(3-2) transition (Fig. Ht) displays 
the be autifully collimated outflow discovered by Bo urke et alJ 
Jl997h . Its PA is ~165°-170°, and its projected collimation fac- 
tor is 4 (as computed on the half-power contour). The cloud 
emission arises between -6 and -3kms _1 . Note that the two 
lobes are well separated spatially, and that there is a small 
amount of red emission in the blue lobe and vice versa. These 
two facts confirm that the inclination of the outflow with re- 
spect to the plane of the sky is close to its semi-o pening an- 
gle (~15°), as already noted by B ourke et alJ Jl997h . The chan- 
nel map also shows hints of a second fainter outflow that is 
a bit tilted relative to the bigger one, with the blue and red 
lobes inverted. This is clearer in Fig.^i, where contours from 
emission for velocities in the [-14,-6] and [-3,5] kms -1 range 
have been plotted. The second outflow is centred on the second 
source of the binary, IRS2. Its PA is between -35° and -30° 
and its projected semi-opening angle ~ 25°. Its projected colli- 
mation factor is 2. No red emission seems associated with the 
blue lobe, and the blue and red lobes are well separated spa- 
tiall y, so the inclination of this outflow is between 25° and 65° 
dCabrit & Bertouil986l) . 

3.2. Temperature enhancement in the lobes 
3.2.1. CO optical depth 

C ls O was observed towards the two bright positions of the 
outflow at offsets (20", -80") and (-40", 140"). We assume 
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isoton ic ratios 12 C/ 13 C = 69+6 and 16 0/ 18 = 557±30 JWilsonl 
\199% . The expected 13 CO/C ls O abundance ratio is then 8.11} |. 
As shown in Fig.^;, the ratio between the antenna temperature 
of the two isotopologues is found to lie in the 8-10 range, so we 
conclude that the 3-2 lines of both isotopologues are optically 
thin. 

From the 13 CO(3-2), we can estimate the optical depth in 
the CO(3-2) transition. We computed the opacity in various ve- 
locity channels (see Tableful. The highest opacity is measured 
on the (0", 0") position. Towards the two lobe positions, the 
opacity is lower: ~ 5 for the (20", -80") offset position and ~ 3 
for (-40", 140"). 

3.2.2. CO excitation 
Table 1. CO LVG modelling. 



Vel channel r 32 N C o T rab , (4 _ 3) /T mb , (3 _2) 

(kms -1 ) (cm -2 ) 



T k 
(K) 



Offset (-40", 140") - Northern lobe (big outflow) 



[-3,0] 
[0,2] 



3.2 
1.4 



1.3xl0 17 
2.9xl0 16 



0.67 
0.65 



>35 
> 20 



Offset (0",0") 



Blue component (IRS1 outflow) 



[-7,-6] 


13 


3.8xl0 16 


0.20 


> 25 




Red component (both outflows) 




[-3,-2] 


<3.5 


<2.4xl0 16 


0.79 


>20' 










>20 2 


[-2,-1] 


<5.5 


<2.5xl0 16 


1.2 


>80 1 










>60 2 


[-1,0] 


<8.1 


<2.5xl0 16 


1.5 


>190> 










>130 2 




Offset (20", -80")- 


Southern lobe 





Blue component (IRS1 outflow) 
[-9,-6] 4.8 2.8X10 17 0.54 > 50 

Red component (emission dominated by IRS2 outflow) 
[-3,0] <3.0 <1.7xl0 16 1.9* > 320 3 

> 300 4 

1.7 < T m b,(4,3)/T ra b,(3,2) <2.1 (lcr). 

1 Assuming N co /Av=3xl0 16 cnr 2 km -1 s. 

2 Assuming N C o/Av=3xl0 15 cm -2 km -1 s. 

3 Assuming N C o/Av=6xl0 15 cm -2 km -1 s. 

4 Assuming N C o/Av=6xl0 14 cm -2 km -1 s. 

Using a standard spherically symmetric large-velocity- 
gradient (LVG) code, we can derive the kinetic temperature in 
the outflow lobes from our CO(3-2), I3 CO(3-2), and CO(4-3) 
observations. The inputs of the code are the CO column density, 
the kinetic temperature, and the density. The CO column den- 
sity can be derived f rom the CO(3-2) optical depth as shown 
bv IChoi et alJ dl993l) . In each velocity channel i, the column 
density is given by: 



Np°(cw" 2 ) 



1.1 x 10' 



T mb ,32Av (Kkms ') 
D(n,T k ) 



T32 



Image Frequency (MHz) 
3. 508 10 5 3.506 1 5 3.504 1 5 3.502 1 5 3.5 1 5 



.5 - 



1.5 — 



» 1 ^ 
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Fig. 2. Methanol spectra observed towards the 3 positions 
shown in Fig. lb, c. The vertical dashed numbered lines in the 
middle panel refer to the following transitions: (1) lij-Orj.o 
[USB], (2) 7 ,7-6o, 6 , (3) 4 , 4 -3_ u [USB], (4) 7_ 1 , 7 -6_ 1 , 6 , (5) 
7o,7-6o,6, (6) 7 +1>6 -6+i,5, (7) 7 2 ,5-6 2 ,4, & 7_ 2 ,6-6-2,5- The inset 
shows the double-peak structure of the 7oj-6o,6 line towards 
the northern position. 



where D(n,Tt) was shown to be equal to 1.5 within a factor of 
2 in the whole range 10 < T k < 200K and 10 4 < n < 10 6 
cm -3 . All the computations are done assuming that the beam- 
filling factor is one for each transition. However, the CO(4-3) 
emission may be less extended than the CO(3-2) emission, so 
the derived 7* and n are lower limits. Table [2 summarises the 
results of this modelling. The constraint on the density is found 
to be not significant (n > 10 3 cm -3 ), and was thus not reported. 

Modelling of the CO emission shows temperature enhance- 
ments from 20 to 50 K towards the IRS1 outflow, which is 
consistent with the temperature derived from methanol anal- 
ysis (see Sect. 3.2.3). It should be noted that both the (0",0") 
and the (20", -80") offset positions comprise emission from the 
small outflow. In particular, the red emission at the (20", -80") 
position is dominated by this second outflow. The derived tem- 
perature is very high, as the CO(4-3) emission is stronger than 
the CO(3-2). The outflow powered by IRS 2 thus appears to 
be warmer than the IRS1 outflow. This could be caused by 
the observation targets including the tips of the IRS2 outflow 
cavity, thereby tracing the present shocks, while toward IRS1, 
the observed position corresponds to a position in the centre of 
the flow, which could trace colder postshock gas. No te that the 
HH321B shock dCorporon & Reipurthlfl997HBourkell200lb is 
encompassed in the beam of our observations at the (20", -80") 
position, and high-resolution IR observations would be very 
useful for determining if this HH object is indeed associated 
with the IRS2 outflow, as suggested by the associa tion of the 
HH320 objects with its northern lobe (Bo urkefcOOll) . 
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Fig. 3. Methanol synthetic spectra. A CH3OH column density of 10 lb cm 2 and a source size of 7" were assumed. 



3.2.3. Methanol statistical equilibrium calculations 

The CH3OH Ik — » 6k band at 338 GHz was observed towards 
the core position and two bright positions of the main out- 
flow with the APEX-2a receiver tuned in the lower side band 
(Fig. |5Jl. This frequency setup allows for simultaneous detec- 
tion in the upper sideband of the 4o4-3_o line at 350.688 GHz 
and of the ln-Ooo transition at 350.905GHz, which have low 
excitation energies (T ow < 20 K) compared to the 7* — > 6a 
band (Ti ow > 44 K). The three positions differ significantly in 
the methanol emission. Towards the central position only weak 
emission is detected from the 4o4-3_i3, the ln-Ooo, and maybe 
the 7()7-6o6 transitions. In contrast, lines from the 7* — > 6k are 
also observed towards the two lobe positions, the lower energy 
transitions toward the northern lobe, and moderately higher en- 
ergy transitions towards the south position. In the (20", -80") 
offset position, the methanol emission is blueshifted, and thus 
associated with the IRS1 outflow. In the (-40", 140") offset po- 
sition, the lines are double-peaked, with a redshifted peak and 
a contribution at the cloud velocity (see inset in Fig. 13, which 
might be associated with emission from the preshock material, 
where methanol is released from the grains, but not accelerated 
yet. 

iLeurini et alJ J2004I) concluded that low-excitation CH3OH 
transitions often trace the density of the gas, while highly 
excited lines are also sensitive to the kinetic temperature. 
Methanol can therefore be used as a tracer of both tempera- 
ture and density. To determine the physical parameters we used 
a spherically symmetric large velocity gradient statistical equi- 
librium code, with the cosmic background as the only radia- 
tion field. Line ratios between different transitions can be used 
to constrain the temperature and density. We assumed that all 
lines are emitted in the same gas; i.e. we suppose that all lines 
have the same beam-filling factor. This assumption is likely to 
be correct, as all the detected transitions have similar lower en- 
ergy (< 75 K, and even between 44 and 75 K, if we only focus 
on the Ik — > 6k band). 

Fig. |3 presents synthetic double-side-band spectra for dif- 
ferent temperature and density conditions, to be compared with 
our observations. We conclude that, towards the (20", -80") po- 
sition, CH3OH comes from a relatively warm (30 < T < 50 K) 



and dense (10 5 < n(H2) < 3xl0 5 cirT 3 ) gas. The temperature 
enhancem ent that we derive i s a bit lower than the range sug- 
gested by iGarav et al.l Jl998l between 50 and 100K, as de- 
duced from analysis of lower energy lines). The transitions they 
used span a narrower range of excitation conditions than the 
presently observed ones, so they are le ss sensitive to tempera- 
ture, as shown bv lLeurini et alJ f2004). The detection of only 
the 7_i7-6_i6 and 7o7-6oe lines towards the red lobe implies 
that the gas temperature must be around 20 to 30 K and not 
denser than a few 10 5 cm 3 . Finally, the non detection of the 
Tk - * 6* lines towards the core position infers a temperature 
~ 10 K (computation for N C H 3 OH=10 I6 cirr 2 ). 2 

4. Conclusion 

From CO and methanol observations, we have constrained 
the density and temperature enhancement in the lobes of 
the highly-collimated outflow of BHR71. Excitation analysis 
points to temperatures between 30 and 50 K, while the den- 
sity is around 10 5 cirT 3 . Our observations have allowed us to 
isolate and determine the structure of the second outflow in 
BHR71 driven by IRS2. From CO excitation analysis, this out- 
flow seems to have very high temperatures. Further observa- 
tions to better constrain the temperature enhancement of this 
second outflow (e.g. the methanol Ik — > 6k band towards the 
blue lobe) would be very valuable. 
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